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in graphene p-n junctions 
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We investigate the transport characteristics of monolayer graphene p-n junctions by means of the 
non-equihbrium Green's function technique. It is shown that thanks to the high interband tunneUng 
of chiral fermions and to a finite bandgap opening when the inversion symmetry of graphene plane 
is broken, a strong negative-differential-conductance behavior with peak-to- valley ratio as large as a 
few tens can be achieved even at room temperature. The dependence of this behavior on the device 
parameters such as the Fermi energy, the barrier height, and the transition length is then discussed. 

PACS numbers: 



I. INTRODUCTION 

Since the discovery of the Esaki tunnel diode [1 , the ef- 
fect of negative-differential-conductance (NDC) has given 
rise to an intense research activity from fundamental as- 
pects of transport to possible applications including oscil- 
lator, frequency multiplier, memory, fast switching, etc. 
|2J . Motivated by the recent development of graphene na- 
noelectronics [MS], this effect has been also investigated 
and discussed in some nanostructures based on mono- 
layer graphene ^ , bilayer graphene [8 , and graphene 
nanoribbons [OtiTT]. However, these studies have shown 
that due to their zero-bandgap, the peak-to- valley ratio 
(PVR) of NDC in 2D-graphene structures is relatively 
small, while, though it can be improved significantly in 
narrow graphene nanoribbons, the results obtained are 
severely affected by increasing the ribbon's width and/or, 
especially, by the edge defects [TOl |TT] . 

A bandgap opening may be a key-step to improve 
the operation of graphene devices. Beyond the tech- 
nique which consists in patterning a graphene sheet into 
nanoribbons , some recent works have suggested that 
a bandgap can open when the inversion symmetry of the 
graphene plane is broken, e. g., by the interaction with 
the substrate [TSHIS], the patterned hydrogen adsorp- 
tion [IB], the adsorption of water molecules [17], or by 
the controlled structural modification [T8^. In particular, 
the experiment reported in [13 has shown that graphene 
epitaxially grown on SiC can exhibit a bandgap up to 260 
meV. Though relatively small compared to that in con- 
ventional semiconductors, it is about ten times greater 
than the thermal energy at room temperature, which has 
stimulated further investigations of graphene-based tran- 
sistors [T9H2T] . In this work, our aim is to look for possi- 
bilities of achieving a strong NDC behavior in monolayer 
graphene p-n junctions. This expectation comes from the 
fact that the appearance of a finite bandgap can result in 



a low valley current, while, the high interband tunneling 
of chiral fermions in graphene may lead to a high current 
peak. Therefore, a large PVR is expected to be achieved 
in these junctions. 



II. MODEL AND CALCULATION 
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FIG. 1: (Color online) Schematic of p-n junction based on 
monolayer graphene (a) and its potential profile along the 
transport direction (b). The n-doped (left) and the p-doped 
(right) regions are used to form a potential barrier in the 
device. The transition between two regions is characterized 
by the length L. 

Graphene has a honeycomb lattice structure with a 
unit cell consisting of two carbon atoms - normally re- 
ferred to as A and B atoms. To describe the charge states 
in the system, a simple nearest neighbor tight binding 
model can be conveniently used [3 , with ac = 0.142 nm 
as carbon-carbon distance, t = 2.7 eV as hopping energy 
between nearest neighbor sites, and = — = A as 
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the onsite energies in the two sublattices. When A = 0, 
this model results in a band structure with zero gap, i.e., 
the conduction and valence bands meet at the K and 
K' points in the Brillouin zone. By making A 7^ 0, the 
inversion symmetry is broken and therefore a finite gap 
opens in graphene's band-structure. The energy disper- 
sion close to the K-point simply writes 



(1) 



where vp = 3act/2h ^ 10^ m/s is the Fermi velocity, k = 
(kx^ky) is the 2D-momentum, and the sign +/- stands 
for the conduction/valence band, respectively. From eq. 
(IT]), the bandgap is determined as Eg = 2 A. To describe 
the excited states around such K-point, one can use the 
following effective (massive Dirac-like) Hamiltonian: 



H = —ihvF {cFxdx + ^ydy) + Acr^ + U 



(2) 



where U is the external potential energy and cr^^y^z are 
the Pauli matrices. The Hamiltonian Q is now used 
to study the transport characteristics of the p-n junc- 
tion schematized in Fig. 1. The p-doped and n-doped 
graphene regions can be generated by electrostatic dop- 
ing [22 , 23 or by chemical doping [24, 25 . The back-gate 
is used to control the carrier density in the graphene layer 
by applying a voltage Vguh- The junction is characterized 
by the potential barrier Uq and the transition length L. 
This length is the size of the region across which the 
charge density changes monotonically from n-type to p- 
type. Though expected to be short [24], it is finite and 
is considered in this work as a parameter. 

In principle, the presence of defects and impurities is 
unavoidable and may result in a substantial amount of 
disorder in the graphene. However, such disorder af- 
fects strongly the electronic characteristics only when 
the graphene sheets are narrow [26]. Thus, we assume 
that it can be negligible in this work where we consider 
the transport in wide systems [27]. Our study addresses 
the ballistic transport through the junction and the non- 
equilibrium Green's function (NEGF) technique (see the 
formulation in ref. [6]) is then used to solve eq. Q and 
to determine the transport quantities. 
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FIG. 2: (Color online) Local density of states (a) and trans- 
mission coefficient (b) in a graphene p-n junction. The tran- 
sition length in (a) L — 20 nm; in (b): L — 5 nm (dotted- 
dashed), 10 nm (solid) and 20 nm (dashed lines). All results 
are for A = 65 meV (A y/E^ + A^ and Ey Tivrky). 
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FIG. 3: (Color online) (a) I — V characteristics of graphene p- 
n junction with different energy bandgaps. Other parameters 
are: Ep = 0.26 eV, Uq = 0.52 eV, and L = 10 nm. (b) and 
(c) illustrate the band diagrams at low bias and in the current 
valley, respectively. 



III. RESULTS AND DISCUSSION 

Before investigating the behavior of the electrical cur- 
rent through these junctions, we plot a map of the local 
density of states and the transmission coefficient as a 
function of energy in Fig. 2(a) and Fig. 2(b), respec- 
tively. The results show clearly three important trans- 
port regions: (i) thermionic when E < Un — A or E > 
[/p + A, (ii) interband tunneling when Un -\- A < E < 
Up — A, and (iii) transmission gap when \E — Ua\ < A, 
where Ua = Un,p denote the potential energies in the 
n-doped and p-doped regions, respectively. The appear- 
ance of transmission gap is essentially due to the fact 



that the longitudinal momentum kx defined from eq. ([T]) 
is imaginary in such energy region, i.e., the carrier states 
are evanescent and therefore decay rapidly when going to 
one of the two junction sides. Note that the same results 
are obtained for A = const even with different Ey and/or 
different A. When A = 0, the transmission gap is just a 
function of Ey and disappears for normal incident parti- 
cles. By rising A, this gap increases and gets its minimum 
value Eg = 2A when Ey = 0. Moreover, due to the ap- 
pearance of evanescent states around the neutral points 
in the transition region. Fig. 2(b) shows clearly that the 
larger the transition length, the lower the interband tun- 
neling. However, it is worth noting that the interband 
tunneling of chiral fermions observed here is still very 
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high in comparison with typical values of the order of 
10~^ observed in conventional Esaki tunnel diodes [28] , 
Therefore, a high current peak and then a large PVR 
are expected to be observed in the considered graphene 
junctions. 

Now, in Fig. 3(a), we display the I — V characteris- 
tics obtained for different energy bandgaps. Note that, 
throughout the work, the current is computed at room 
temperature. It is shown that the NDC behavior ap- 
pears more clearly and its PVR increases when increas- 
ing the bandgap though the current peak reduces. The 
form of the I — V curve can be explained well by looking 
at the diagrams in Fig. 3(b,c): at low bias (3(b)), the 
contribution of interband tunneling processes makes the 
current rising; when further rising the bias (3(c)), the in- 
terband tunneling is suppressed due to the transmission 
gaps and therefore the current decreases; when the bias 
is high enough, the contribution of charge carriers in the 
thermionic region makes the current rising very rapidly. 
Because the transmission gap increases, the current, es- 
pecially in the valley region, decreases when increasing 
the bandgap. This results in a large PVR of NDC, for 
instance, about 46 in the case of Eg = 260 meV shown 
in Fig. 3(a). 
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FIG. 4: (Color online) I — V characteristics of graphene p-n 
junction with different Fermi energies (a) and different po- 
tential barriers (b). Other parameters are: Eg = 260 meV, 
Ef = Un + Uo/2 in (b), Uo = 0.52 eV in (a), and L = 10 nm, 
where Un is the potential energy in the n-doped graphene. 
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FIG. 5: (Color online) I — V characteristics of graphene p-n 
junction with different transition lengths. The inset shows 
the evolution of peak-to- valley ratio with respect to L. Other 
structure parameters are: Eg — 260 meV, Ef — 0.26 eV, and 
Uq = 0.52 eV. 

To evaluate the possibilities of having a strong NDC 
behavior in the junctions, we consider the behavior of the 
I — V characteristics with different Fermi energies in Fig. 
4(a) and different barrier heights in Fig. 4(b). Fig. 4(a) 
shows that the current valley shifts to the low bias region 
while the current peak decreases when Ep goes from the 
value of Um = Un to the top of the valence band 

in the p-doped region or to the bottom of conduction 
band in the n-doped one. Simultaneously, the current 
is a symmetric function of Ep around Um- This result 
leads to the fact that the strongest NDC is achieved when 
Ep = Um and the rectification behavior around the zero 
bias can be observed when Ep goes away from such value, 
e.g., see the cases of Ep — 0.13 or 0.39 eV in Fig. 4(a). 
The latter behavior is essentially due to the role of the 
transmission gap (as illustrated in the diagrams in Figs. 
3(b,c)), which induces a strong reduction of current in the 
low-positive bias region. Besides, when decreasing the 
barrier height, we find that the transmission gap around 
Up moves downward in energy, and as a consequence, the 
interband tunneling and the current peak are reduced. 
Indeed, this point is clearly illustrated in Fig. 4(b), where 
Ep is chosen to be U^ -^Uq/2 to achieve the strongest 
NDC. Moreover, it is shown that the current valley shifts 
to the low bias and its value is reduced, which finally 
results in an increase of the PVR when decreasing Uq. 
For instance, the PVR is - 30, 46, and 79 for /7o = 0.6, 
0.52, and 0.44 eV, respectively. To obtain the best device 
operation (either a high current peak or a large PVR) , the 
study suggests that the barrier height about two times 
greater than the energy bandgap should be used. 

Finally, the I — V characteristics in the junctions with 
different transition lengths are displayed in Fig. 5. It is 
shown that because of the suppression of interband tun- 
neling (as in Fig. 2(b)), the current peak is reduced when 
increasing L. Therefore, as seen in the inset of Fig. 5, 
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the PVR decays as a function of the length L. However, 
the PVR as large as about 10 is still obtained even for L 
up to 20 nm. This suggests the possibility of achieving 
a very large PVR with a short transition length, which 
can be realized by controlling the device geometry, e.g., 
by appropriately reducing the gate dielectric thickness in 
the case of the electrostatic doping [23^, or by using the 
chemical doping to generate the junction as mentioned 
in ref. |24]. As expected above the large PVR, e.g., ~ 
123 obtained for L = 5 nm in this work, is much higher 
than for conventional Esaki tunnel diodes where its high- 
est reported value is just about 16 (see in ref. [29] and 
references therein). 

IV. CONCLUSION 

Using the NEGF technique, we have investigated the 
transport characteristics of monolayer graphene p-n junc- 
tions. Even at room temperature a negative-differential- 
conductance with peak-to-valley ratio as large as a few 
tens can be observed in these junctions thanks to a fi- 



nite bandgap and to the high interband tunneling of 
chiral fermions. The dependence of this behavior on 
the device parameters was analyzed. It is shown that 
the strong negative-differential-conductance behavior is 
achieved when the Fermi level is in the center of the 
potential barrier, the barrier height is about two times 
greater than the bandgap, and the transition length 
is short. We hope that obtained results can be help- 
ful for designing efficient room-temperature negative- 
differential-conductance devices based on finite-bandgap 
graphene. 
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